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Про визначення контакту у модельному класі Сретенського 


(Представлено академіком НАН України В.І. Старостенком) 


Сформульована обернена контактна задача гравіметрії для контактної поверхні в класі Сретенського. 
Доведено розривність вертикальної похідної сили тяжіння, яка фігурує в її правій частині. Ця особли- 
вість використана для редукції задачі до альтернативної постановки у вигляді інтегрального рівняння 1- 
го роду. Шляхом різницевого аналізу меж відповідних інтегралів виведено ефективний спосіб послідовно- 
го обчислення контактів. Задача узагальнюється на п меж і зберігає коректність в межах класу Срете- 
нського. 


Вступ. З початку ХХІ ст. відбувається неперервна зміна основної парадигми теорії інтерпрета- 
ції потенціальних полів |1|. Вона супроводжуєтеся переглядом ключових методологічних (2), апа- 
ратно-програмних (|3| та чисельних |4| засад тлумачення даних цих полів. В рамках творення чисе- 
льного базису цієї нової парадигми особливої ваги набуває розробка адекватного геофізичній пра- 
ктиці апарату математичного моделювання геофізичних подів, адаптованого до обробки масивів 
даних великої розмірності, у тому числі - на комплексах паралельних обчислень. 

Ключове місце в рамках нової методології належить аналітичним апроксимаціям середовища й 
поля. Теоретичне підгрунтя такого підходу викладено в роботі Страхова (2). В рамках нової мето- 
дології активно розробляються конструктивні аналітичні апроксимації геологічного середовища. 
Зокрема, введені в науковий обіг апроксимації на основі типових блоково-шаруватих або сферич- 
но-циліндричних елементів |5|. Постульовано нагальну потребу доповнити ці конструкції апрок- 
симаціями потенціальних полів |6|. Це, певною мірою, було зроблено у праці |7|, в якій запропо- 
новано новий математичний базис (диференціальне рівняння) для опису аномалій сили тяжіння. Та- 
кож вивчаються апроксимації рельєфу поверхні спостережень шляхом параметризації заданими 
гладкими функціями |З. 

Продовжуючи актуальні розробки аналітичних апроксимацій середовища, ми пропонуємо нову 
апроксимаційну схему для визначення контакту двох однорідних тяжіючих тіл, які належать до 
класу Сретенського, за заданим розподілом вертикальної компоненти аномалій сили тяжіння. Зау- 
важимо, що аналітичні конструкції методу дозволяють реалізувати розпаралелювання обчислень 
за рахунок рекурентності виразів для розв'язання відповідної прямої задачі. 

Відзначимо, що постановка і розв'язання прямих та обернених задач граві- і магнітометрії у 
постановці для класу тіл Л.М. Сретенського започатковані в працях Булаха |9). Ці постановки ви- 
писані для чисельної моделі як із одного |10), так і кількох (11| рудних тіл. В силу вдалої та зруч- 
ної параметризації вони зручні для практичного застосування в комплексі інтерпретації даних руд- 
ної геофізики. Однак формулювання, власне, самого класу Сретенського моделей аномальних тіл 
дозволяє виписати пряму задачу гравіметрії у структурній постановці, чого досі не було зроблено. 
Поширення розв'язку цієї задачі на кілька контактних поверхонь відповідатиме сучасній практиці 
моделювання складного геологічного середовища за комплексом геофізичних даних. 

Постановка задачі. Нехай на деякій горизонтальній площині 2 - 0 (вісь 2 спрямована вниз) 
задані значення вертикальної похідної (/, (2.0) потенціалу сили тяжіння, зумовленого "нескінчен- 


ним" шаром С (видовжується в обидва боки за профіль спостережень), витягнутим уздовж осі у. 
Нехай цей шар складається із двох пластів із постійною густиною с « сопзі, розділених горизонта- 
льною контактною поверхнею Н. Цей об'єкт належить класу двовимірних тіл, для яких існує сере- 


дня площина 2. (вона проходить через тіло так, що будь-який перпендикуляр до неї перетинає по- 
верхню тіла лише у двох точках, по різні сторони від площини). Контактна поверхня Н ундулює 


навколо середньої площини Ук , не надто ухиляючись від неї, У смислі метрики бананового прос- 


тору Н п Хв - паї М В частинному випадку Н може співпадати З ДИ й Назвемо цей модельний 
хер 


клас тіл - класом Сретенського 5г(, Д), (структурний варіанту), де 1 означає сталу густину, а Р - 


замкнену область (поверхню), виповнену тяжіючими масами. 
Зауваження 1. Якщо площина спостережень (денна поверхня) не є горизонтальною, її можна 


апроксимувати будь-якою конструкцією 27 з н(», у), у тому числі, однією із пропонованих в |З). 


Зауваження 2. Визначення класу 5г(, р) через поняття середня площина не є конструктивним, 


оскільки застерігає лише від вжитку "глечикоподібних?" тіл і ніяк не регламентує спосіб визначен- 

ня цієї середньої площини, наприклад, для множини багатокутників, серед яких можуть трапитись 
і такі тіла, як зображене на рис. І. 

На нашу думку, подолати амбівалентність означення класу 


5, р) можна, залучивши диференціальні характеристики конк- 


ретного сімейства тіл, як, наприклад, зроблено в (8). Втім, подіб- 
ні міркування виходять за рамки статті. Модель горизонтального 
циліндра апріорі задовольняє вимоги класу. 


Уточнимо характеристику класу 51, Б). Сполучаючи із се- 
редньою площиною аномального тіла координатну площину 
201 , аналітично опишемо область І), зайняту тілом, так: 


раКа бака єв, С(б)« с «СЕ, а) 
Рис. 1. Модель рудного тіла, для Де 90406 (2) «со - певні контактні поверхні. При цьому ме- 


якої середня площина у смислі Л. жу ОР області )) можна представити як об'єднання двох конту- 
Сретенського є неоднозначною. й Й 
рів 2, 00Ю, у вигляді: 


др, «(с ,б)гач є ССС) і 212). (2) 
Розглянемо підклас 5, 0.) класу 5, р), коли а - -00, Б - 00, аосі координат 2 іх, 0 і 
2 паралельні. Тоді, не порушуючи загальності міркувань, можна представити область 1), і її межі 


у аналітичному вигляді: 


ра Ще є): -оса«огт)«/х« 12) 
др г а, с): -о«б со, 19 92), із 12) (3) 
Нехай підклас 5, г) тяжіючих тіл складається із двох контактних поверхонь 
ба) ока «сі зі2, 4 
розділених середньою площиною 2 - 2,20. 


Поставимо таку задачу: визначити контакти (4) за значеннями вертикальної похідної (/, (2.0) 


потенціалу сили тяжіння. 
Аналітичні властивості аномалій. Розв'язати цю задачу допоможе аналіз властивостей похід- 


них потенціалу сили тяжіння (/ (х, Уа ї; Представимо ці похідні у вигляді 
90 (Сх а) де 
а2о азас, 
зе а ря 


а ие 2-0 АС, (5) 


де с (6-х) (с -2) 
де, в даному випадку, 


рес): ооо, 00) (6) 


Внаслідок подання (6), подвійні інтеграли в формулах (5) можна переписати у такому вигляді: 


о С(5) 
22) зла Їає | кі? й ді 


ох ЦИ Креру зро 
бпбоз) зо ЛО . 
вн ср зраз 


Зважаючи на результати, отримані в праці |(12|, звідси, після ряду нескладних аналітичних 
перетворень, будемо мати 





90, 1) 
Ох 


сах 


цого ша ар 


2д ас (х) «0 
о 2 2 
б «). То Пп е -хя я «Ї 42 44 2діо|С(х)- 2 оч хг. 
2 в 
-24/аг (х), С(х)є 2 
Отже, функція (/, (ж, 1) щодо 27 є неперервною, в той час як функція (/ йРо 1) щодо 2 - розривна. 
На підставі дослідження фундаментальних властивостей (7) контактів (4), за методикою аналізу 


гармонічних функцій, розвиненою в (12), одержуємо при 2 -» - 0 такі вирази зовнішніх похідних 
гравітаційного потенціалу 





во (з 0) аа Їжа Я Чо, 


і-і -оо 


б (2-х |С9(х 
Лимонад о 


У той же час на середній площині 2 - 2, ( С) «кр Су), -о00« ХХ «00) з тих же співвідно- 





шень (7) знаходимо, що 
90, за) 
дх 
дих, зо) 
дг 
«Уа; і і | ЧЕО)-а| аб . 0) 
іні о -Х 


- (2 Ї 
- І (2)- Хо 
Спосіб розв'язання. Очевидно, що значення потенціалу на поверхні (8) і на середній площині 
(9) розбігаються. Це нескладно було передбачити, виходячи з виразу (7). Однак, у цьому протиріч- 
чі закладена цікава можливість для розв'язання поставленої задачі. 
Дійсно, знайдемо внутрішню межу похідної, обчисливши інтеграл Пуассона 
2050) Я і дп о) Ко ан (10) 
Ли 2 2 ; 
Ог л 07. (п - х) а 
При цьому представимо підінтегральну функцію у наступному вигляді 
ди с 25 У ї 
-21о, Пра у асас з 
п 


з Час 


(9 
зе із час, 


щ дос (2) -Ол/о, РАЗТРУ ТА ан 


аз Ра 
2-х) (а -2) 
для тяжіючої області С; з Хе, б. у: -ох б «оба ко Фе Обчислюючи інтеграл (10), оде- 
ржимо, з урахуванням наведеного вище зображення, такий вираз: 
о, (2.0) (1 (0) 
ба т тек атома (х)- 
(2-х) «ЕОООЇ 
/ 2 най. (п) 
їаЇ 2-х) РОЇ 
Враховуючи значення зовнішньої межі (8), та величину б (11), знайдемо різницю меж: 
УкАТО 
ги) 8) она їжі у «Кеооо, 
ді дз Ра У 2 
Щоб спростити подальші викладки, позначимо цю різницю потенціалів через функцію Ух), 
1 (оо. (с 0) д0.(х,0) 
4ліо, Ох 02 


аа 











де У(х)- 


. Відтак, ми отримали нелінійне інтегральне рівняння 


1 25 Р х) С (х) Її 
с9(х)- б ас -У(х). (2) 
2л їі уче (2) Її 
для визначення першого наближення контакту 20 С), - оо «хор, Розв'язання подібних 


рівнянь започатковане у вже заданій вище праці |12). 
Р я 4 1 
Визначивши із останнього рівняння функцію 0 ( 1), зможемо обчислити наближення поля 


задо - 2доуб"(х)- є Їи й - 5 . вна 46. ал) 


А це, у свою чергу, дозволяє знайти в "чистому вигляді" ефект від наступного контакту у вигляді 
різниці граничного та попереднього значень вертикальної похідної потенціалу: 


20, (о0) д0,(50) д0 (0) 
да да ді 


Маючи цю величину, можна з рівняння, аналогічного рівнянню (13.1), тобто 


5 7. (Є-ху з СО) 90. (х0) 
2дісьс! з го, Па б ра о оз (13.2) 


2 2. - 
обчислити чергове наближення контакту 2 ) ша с Ох), -о«х«со0. Таким чином, ми знайшли 





"верхню" 1 "нижню"? межі області (1), що й розв'язує поставлену задачу (4). До речі, власне, самі 
контакти С Фе ) можна апроксимувати, серед іншого, конструкціями типу наведених в праці (111. 

Зауваження 3. Наведені вище математичні викладки легко узагальнюються на випадок п меж 
10)- С ж), -о«х «00, якщо між кожною парою цих меж можна провести середню площину 


(тобто якщо ці межі не виходять за межі класу 5, р)). 


Висновок. Введено в науковий обіг нові апроксимаційні конструкції (12) й (13.1-2), за допомо- 
гою яких можна ефективно відновлювати будову шаруватого геологічного середовища, зображе- 
ного парами (верхня і нижня межі) контактів, обмежених властивостями класу Сретенського. 


Доведено, що на класі 5, р) послідовні наближення 0 2) однозначні й стійкі (збігаються 
0 : - 
до заданого точного розв'язку 0 | 2), якщо похибка вхідних даних не перевищує 390 від макси- 
мальної амплітуди аномалій Ї/, (х,0) Г131. 

Дискретизацію інтегральних рівнянь (12) й (13.1-2) доцільно здійснити за допомогою відомих 
методів скінченновимірної редукції інтегральних рівнянь 1-го роду |14|. Алгоритм розв'язання за- 
дачі пройшов первинну апробацію на працездатність. У перспективі передбачається поширення 
аналітичних викладок на тривимірний випадок та апробація алгоритму на практичних матеріалах. 
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Ю.І. Дубовенко 

ПРО ВИЗНАЧЕННЯ КОНТАКТУ У МОДЕЛЬНОМУ КЛАСІ СРЕТЕНСЬКОГО 

Сформульована обернена контактна задача гравіметрії для контактної поверхні в класі Сретенського. 
Доведено розривність вертикальної похідної сили тяжіння, яка фігурує в її правій частині. Ця особливість 
використана для редукції задачі до альтернативної постановки у вигляді інтегрального рівняння 1-го роду. 
Шляхом різницевого аналізу меж відповідних інтегралів виведено ефективний спосіб послідовного обчисле- 
ння контактів. Задача узагальнюється на п меж і зберігає коректність в межах класу Сретенського. 

Ключові слова: гравіметрія; обернена задача; аналітична апроксимація; клас Сретенського; відновлення 
контакту; контактна поверхня; інтегральне рівняння; гравітаційний потенціал; математична модель; середня 
площина. 


Ю.Й. Дубовенко 

ОБ ОПРЕДЕЛЕНИЙ КОНТАКТА В МОДЕЛЬНОМ КЛАССЕ СРЕТЕНСКОГО 

Сформулирована обратная контактная задача гравиметриий для контактной поверхности в клас- 
се Сретенского. Доказана разрьвность вертикальной производной силь тяжести, стоящей веб пра- 
вой части. Зта особенность использована для сведения задачи к альтернативной постановке в виде 
интегрального уравнения 1-го рода. Путем разностного анализа границ соответствующих интегра- 
Лов вьіведен зффективньюй способ последовательного вьтчисления контактов. Задача обобщаєтся 
на п границ и сохраняєт корректность в рамках класса Сретенского. 

Ключевь(е слова: гравиметрия; обратная задача; аналитическая аппроксимация; класс Сретен- 
ского; контактная поверхность; интегральное уравнениєе; гравитационньм8й потенциал; математиче- 
ская модель; средняя плоскость. 
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ОХ ТНЕ рЕМ5ІТУ ІХТЕВКАСЕ РЕКІХІТІОКХ У(ІТНІХ 5КВЕТЕ5КП МОРЕІЇ, СІ.А55 

ТРе сопіасі ргобіет ої огамісу іпуег5іоп Гог деп5ісу іпіегіасе ууїріп Ше Згегеп5Кії сіа55 15 розі(шіаєеа. 
А дїзсопійпиїу ої Де єгаміїу уегіїса! дегіуайуе, убісП арреаг8 їп Ше ргобіепія гід Бапа, 15 ргоуеа. 
ТБіз Тваиге 15 цзед їо гедисе їбе ргобіет їо аПегпайуєе 5(акетепі іп їп пе Їогт ої Бе І" Кіпа іпіевтаї! 
едиайоп. Ву пеапз ої гезідца! апаіузіз ої (пе Боппаагіе5 Гог пе соггезропаіпя іпіезтаїз, ап ебісіепі (е- 
српідше Гог 5иссе8зіуе саїсціаноп ої Бе деп5іїу іпіепіасе 15 дедисей. ТПе ргобіет 15 сепегай7ед опіо п 
Боипаагіез са5е, ме рге5егутпя, ї(5 соггесіпе85 м'іііп пе Згекеп5Кії сІа55. 

Кеу ууогая: дгауітеїгу; іпуегвіоп; апаїуйсаї! арргохітайоп; 5гегеп5Кії сіа55; Фепзіу іп(іегіасе; іпіеє- 
га! едиайоп; єзгаміїу рогепца!; плаїретайсаї птоде!; Ше теап 5игіасе. 
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